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ABSTRACT 

We obtained spectra of 60 red, starlike objects (E < 18.8) identified with FIRST 
radio sources, Siaghz > 1 mJy. Eight are QSOs with redshift z > 3.6. Combined 
with our pilot search (Benn et al 2002), our sample of 121 candidates yields a total of 
18 z > 3.6 QSOs (10 of these with z > 4.0). 8% of candidates with S 1A gh z < 10 mJy, 
and 37% of candidates with Siaghz > 10 mJy are QSOs with z > 3.6. The surface 
density of E < 18.8, Siaghz > ImJy, z > 4 QSOs is 0.003 deg -2 . This is currently the 
only well-defined sample of radio-loud QSOs at z m 4 selected independently of radio 
spectral index. The QSOs are highly luminous in the optical (8 have Mg < -28, qo = 
0.5, Hq = 50 kms _1 Mpc _1 ). The SEDs are as varied as those seen in optical searches 
for high-rcdshift QSOs, but the fraction of objects with weak (strongly self-absorbed) 
Lya emission is marginally higher (3 out of 18) than for high-rcdshift QSOs from 
SDSS (5 out of 96). 

Key words: quasars: general - quasars: emission lines - radio continuum: galaxies - 
early Universe 



1 INTRODUCTION 

The evolution with redshift of the space density of QSOs 
places strong constraints on the abundance of massive ob- 
jects at the earliest cosmological epochs. Higti-redshift QSOs 
can be found efficiently by colour selection from very large 
samples (e.g. the Sloan Digital Sky Survey, Anderson et al 
2001, see Benn et al 2002, hereafter Paper I, for a summary 
of previous searches), but space-density measurements based 
on such selection could be biased by e.g. redshift-dependent 
dust extinction. 

Any such bias is reduced when selecting in the radio, 
but searches for tiigh-redshift radio QSOs have so far concen- 
trated on radio-bright objects with flat radio spectra, which 
yields small samples with relatively low surface density on 
the sky. E.g. Snellen et al (2001) and Hook et al (2002) 
sought red starlike optical counterparts of S$ghz > 30 mJy 
flat-spectrum radio sources, and identified a total of 8 z > 4 
QSOs. Here, we aim to identify high-redshift radio quasars 
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without any bias in radio spectral index, and with higher 
surface density on the sky, through spectroscopy of red star- 
like counterparts of FIRST radio sources, with Siaghz > 1 
mJy. Combined with our pilot search (Paper I), this yields 
18 z > 3.6 radio-loud QSOs (9 previously known, mainly 
from multi-colour optical searches). This includes the largest 
sample to date of radio-selected QSOs at z > 4. 



2 SAMPLE 

The selection procedure is a refinement of that described 
in Paper I. We sought red, starlike optical identifications of 
the 722354 sources in the July 2000 edition of the FIRST 
catalogue of radio sources (7988 deg 2 , 7 h < RA < 17\ -5° 
< dec < 57°, White et al 1997, Siagb* (peak) > 1 mJy), 
satisfying the following criteria: 



1 All 18 QSOs are radio-loud according to the criterion of Gregg 
et al 1996, i.e. P 1A GHz > 25.5; see Table 1). 
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Table 2. Candidates which are not z > 3.6 QSOs 
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Stocke et al (1991), BL Lac, z = 0.346 

1.893 0.010 CIII 1909, CIII 2326, Mgll 2798 A 
Boyle ct al (1990), BAL, z = 2.18 
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0.1906 0.0005 OIII 4959/5007, H a 

Extended on JKT image, no spectrum obtained 

Extended on JKT image, no spectrum obtained 



Columns 1-6 and 8-10 are as in Table 1. Column 7 indicates spectral type (see caption of Fig. 3). For 3 of the QSOs, rcdshifts are 
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(i) starlike optical counterpart in the APM (Irwin et al 
1994) catalogue of the POSS-I survey (which includes 7030 
deg 2 of the FIRST survey), lying within 1.5 arcsec of the 
FIRST radio source; 

(ii) POSS-I/APM E < 18.8; 

(iii) POSS-I/APM O - E > 3, or no measured O mag 
(plate limit O « 21.7); 

(iv) starlike on POSS-II images (this excludes 80% of can- 
didates satisfying the above criteria); 

(v) red also in the Minnesota APS (Pennington et al 
1993) catalogue of the POSS-I plates, i.e. not omitted from 
the APM catalogue because of confusion on the blue plate, 
see Paper I (this excludes 16% of the candidates satisfying 
the above criteria). 

This yielded 294 candidates, of which we visually classified 
194 as 'starlike' on POSS-II and 100 as 'possibly starlike'. 

CCD images were obtained in good seeing of about one 
third of these, using the Loiano 1.5-m telescope of the As- 
tronomical Observatory of Bologna. On these images, 25 of 
the 70 'starlike' on POSS-II, and 28 of the 33 'possibly star- 
like' candidates were extended. We therefore excluded from 
our sample the 100 'possibly starlike' objects, implying an 
incompleteness of 100 x 5/33 = 15.2 stellar candidates. 

The final list of 194 candidates is therefore expected to 
include 45/70 x 194 = 124.7 truly starlike images, out of 
a total of 124.7 + 15.2 = 139.9 amongst the original list of 
294 candidates. Therefore, it includes 89% (= 124.7 / 139.9) 
of the red QSO identifications of sources in this edition of 
the FIRST catalogue. The selection criteria differ from the 
preliminary criteria used in Paper I only in going slightly 
deeper, i.e. E < 18.8 rather than E < 18.6. 



3 OBSERVATIONS AND REDUCTION 
3.1 Optical Spectroscopy 

Observations of 54 of the 194 candidates were reported in 
Paper I. The 55th object reported there, 1349+38, is ac- 
tually detected on the Minnesota APS scan of the POSS-I 
blue plate (i.e. it is not red), and thus does not meet selec- 
tion criterion (v) above. Spectra were obtained of 60 more 
candidates using the IDS spectrograph on the Isaac Newton 
Telescope during near-photometric nights (occasional light 
cirrus) in 2001 May. The spectrograph was used with the 
R150V grating, yielding spectra with dispersion 6.5 A/pixel, 
usually centred at 6500 A. For details of the observing and 
data-reduction, refer to Paper I. The results are summarised 
in Tables 1 [z > 3.6 QSOs) and 2 (other objects). Table 2 
includes an additional 3 QSOs which were not observed, but 
whose redshifts were obtained from the literature. Spectra 
of the QSOs with redshift z > 3.6 are shown in Figure 1. 

To test the completeness of the search for z > 4 QSOs, 
we observed an additional 28 candidates with 2.0 < O — E < 
3.0. None had z > 3.6. The spectra of two objects with z = 
3.2, 3.4 are given in Figure 2 (details included in Table 1). 

In total, we have spectroscopic information for 117 of 
the 194 candidates with O — E > 3 (54 observed Paper I, 
60 observed here, 3 redshifts z < 3.6 from the literature). 
In addition, we were able to classify 4 objects as galaxies on 
the basis of JKT imaging (see below). The colour-magnitude 
distribution of these 121 candidates is shown in Figure 3. 



The distribution of our 18 z > 3.6 QSOs in radio flux 
density and E mag is compared with that from other radio 
searches for high-redshift QSOs in Fig. 4. In the optical, 
these QSOs are highly luminous, most with Mab (1450A) < 
-27. 

Of the 121 observed candidates, 18 (15%) are QSOs 
with z > 3.6, and 10 of these have z > 4. 

3.2 Optical imaging 

Some of the 194 candidates were imaged through a Harris R 
filter with the 1.0-m Jacobus Kapteyn Telescope (JKT) in 
service time, in order to identify extended objects. Spectra 
were later obtained of most of these candidates at the INT 
(above), but the 'C (galaxy) classifications of 4 objects in 
Table 2 are based on JKT imaging alone. 

In addition, Harris R— and J— band photometry was ob- 
tained with the JKT (see Table 1) for 9 of the z > 3.6 QSOs, 
primarily to allow an independent check of the derivation of 
Mab (1450A) absolute magnitudes (see Vigotti et al 2003). 
Conditions were photometric, except for light cirrus during 
observation of 1639+43. Otherwise, the estimated errors on 
the R, I mags are < 0.05 mag. For smooth continuum spec- 
tra, one would expect POSS-I E and JKT R mags to differ 
by < 0.1 mag (Humphreys et al 1991), at a given epoch. The 
observed R — E values range -0.8 to 1.4 (median 0.7), mainly 
because the R filter covers a wavelength range a factor of 

3 larger than does the E filter, and includes more of the 
depressed continuum blueward of Lya. The range of R — E 
values may also reflect changes in QSO luminosity during 
the 50 years between the POSS-I and JKT observations. 

3.3 Radio observations 

The 13 z > 3.6 QSOs with Siaghz > 3 mjy were observed 
at 4.85 and 10 GHz, with the Effelsberg radio telescope, dur- 
ing the night of 2002 Jan 13. The sources are all unresolved 
by the telescope beam, FWHM 143 arcsec and 69 arcsec at 
these two frequencies. The sources were observed by scan- 
ning the telescope in right ascension and declination, with 
scan speeds 30 arcmin and 10 arcmin per minute respec- 
tively, scan lengths 15 arcmin and 7 arcmin, with the total 
number of scans adjusted to the total flux density of each 
source. The flux densities (scale of Baars et al 1977) and 
spectral indices (between 1.4 and 4.85 GHz) are given in 
Table 1. The median spectral index is a = -0.3 (S v oc v a ). 

4 RESULTS 

4.1 High-redshift QSOs 

The sample comprises 18 z > 3.6 QSOs (Table 1). Below we 
derive the surface density on the sky, and we comment on 
the spectra of individual objects. 

4-1.1 Surface density 

The fraction of high-redshift QSOs amongst the candidates 
does not depend on optical magnitude, but depends strongly 
on radio flux density (Table 3): 8% for Siaghz < 10 mjy, 
37% for Siaghz > lOmJy. 73 of the 194 candidates have not 
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Figure 1. Spectra of the 8 z > 3.6 QSOs amongst the 60 newly-observed candidates (see Paper I for those previously discovered). 
Spectral features are labelled at wavelengths corresponding to the quoted redshift, assuming rest-frame wavelengths in A of 1216 (Lya), 
1240 (NV), 1302 (Ol/Sill blend), 1400 (SilV/OIV] blend), 1549 (CIV). The spectra have not been corrected for terrestrial atmospheric 
absorption, notably at 7594 and 6867 A (A and B bands). The ticks below the spectrum of 0918+06 indicate a z = 4.140 absorption 
system (1260 A Sill, 1302 A Ol/Sill, see also Snellen et al 2001), i.e. velocity 2330 km s" 1 relative to the QSO. Those below the spectra 
of 1303+00 and 1412-01 indicate absorption at z = 3.600 and 3.624 respectively, i.e. outflow velocities 2600 and 6160 km s _1 . 
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Figure 2. Spectra of 2 z > 3 QSOs discovered when observing a sample of bluer, 2.0 < O — E < 3.0, candidates (i.e. not included in 
the main sample of 121 observed candidates). The ticks below the spectrum of 1624+37 indicate an absorption system at z = 3.327, i.e. 
velocity 2970 km s" 1 relative to the QSO. This QSO also exhibits a CIV broad absorption line (BAL) with outflow velocity up to 29000 
km s _1 . 



Table 3. Fractions of candidates which are QSOs with z > 3.6 
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Columns are: (1) FIRST 1.4-GHz integrated flux density (mjy), (2-4) number of z > 3.6 QSOs / total number of candidates, for each 
range of POSS-I E magnitude, (5) sum over all ranges of E, and (6) sum expressed as a percentage. 



been observed, 67 with Siaghz < 10 mjy, 6 with Siaghz > 
10 mjy. On the basis of the above high-redshift-QSO frac- 
tions, we expect to have missed 5.95 z > 3.6 QSOs in 
the sample of 194 candidates, i.e. the completeness of our 
search for candidates satisfying our selection criteria is (18 
/ (18+5.95)) = 75% over the 7030 deg 2 overlap between 
FIRST and the APM/POSS-I catalogues, within 7 h < RA 
< 17 . Multiplying by the 89% completeness due to filtering 
out 'possibly-starlike' candidates (Section 2), our search is 
67% complete, i.e. the surface density of z > 4 QSOs (10 
QSOs E < 18.8, Siaghz > 1 mjy) found with this tech- 
nique is 0.0021 deg -2 . 

This is about twice that of recent searches for high- 
redshift QSOs amongst the red starlike optical counterparts 
of (radio-brighter) flat-spectrum sources in the northern 
(Snellen et al 2001; 4 z > 4, E < 19, 0.0006 deg" 2 ) and 
southern (Hook et al 2002; 4 z > 4, S 5 ghz > 72 mjy, R < 
22.5, 0.001 deg -2 ) hemispheres. The surface densities mea- 
sured for the FIRST sources are higher in part because of 
the much fainter radio flux-density limit, and in part because 
steep-spectrum sources are not excluded. 

In the sample of Hook et al (2002), the ratio of the 
numbers of E < 21 and E < 19 z > 4 QSOs is 7 / 4 = 1.8. 
Combining this with our above result, we predict the surface 
density of E < 21, Siaghz > 1 mjy QSOs with z > 4 to be 
0.004 QSOs deg -2 . 

The surface- density calculations above account for in- 
completeness due to rejection of candidates not considered 
definitely starlike on POSS-II, and due to only 121 out of 194 



candidates being observed. In calculating the overall com- 
pleteness of a search for z > 4 QSOs with Siaghz > 1 nijy, 
one must also take into account: 

(i) the completeness of the FIRST catalogue at low flux 
densities, 83% for the flux-density distribution of Table 1 
(Prandoni et al 2001); 

(ii) the completeness of the APM catalogue of POSS-I, 
84% (see Vigotti et al 2003, section 2.3); 

(iii) that the radio-optical separations of some QSOs may 
exceed 1.5 arcsec due to measurement errors, completeness 
99%, given combined radio-optical rms 0.5 arscec; 

(iv) that some QSOs may be missed because they do not 
coincide with peaks of radio emission (e.g. near the mid- 
points of double radio sources), completeness 98% (see Vig- 
otti et al 2003, section 2.3); 

(v) that some z > 4 QSOs may have O — E < 3, com- 
pleteness ~ 100% (Section 3.1, see also Vigotti et al 2003, 
figure 1). 

The completeness due to these five factors combined is 0.68, 
i.e. the corrected surface densities of z > 4 QSOs with 
Siaghz > 1 mjy is 0.0031 deg" 2 for E < 18.8, and is 0.006 
deg" 2 for E < 21. 

Since ~ 10% of QSOs are radio-loud (e.g. Sramek & 
Weedman 1980; and Petric et al 2003 for consistency at z ~ 
5), we expect the density of radio-quiet + radio loud z > 4 
QSOs with E < 21 to be ~ 0.06 deg" 2 . E.g. Kennefick et 
al (1995a,b) find 0.03 deg" 2 with r < 20. The SDSS search 
(e.g. Anderson et al 2001) reaches a higher surface density, 
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Figure 3. (a) Distribution in colour and magnitude of the 194 starlike candidates, (b) Expanded view of the lower-right corner of Fig. 
3a. The symbols for spectral type (121 objects) are: number = redshift of QSO with z > 3.6, 'Q' = QSO with z < 3.6, 'S' = emission-line 
galaxy (starburst or Scyfert), 'G' = radio galaxy, '?' = probable radio galaxy (definitely not high-rcdshift QSO, but might be a star), 
'*' = star. Large font indicates S\,&ghz > 10 m Jy; small font S\,iGHz < 10 m Jy- Objects with no O magnitude have been plotted at 
the nominal O plate limit (dashed line) with O — E = 21.7 - E. The unobserved 73 = 194 - 121 objects are plotted as circles. In Figure 
3a, the positions of a few points have been adjusted to reduce overlap of labels. In Figure 3b, the E mags have been randomised ± 0.05 
mag for the same reason. 
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~ 0.1 deg 2 because the i and a bands are also used, i.e. no 
effective limit in r band. 



4-1-2 Lya emission 

The rest-frame Lya+NV equivalent widths of the 18 z > 
3.6 QSOs are given in Table 1. The values are approximate, 
rms ~ 30%, due to the difficulty of estimating the slope of 
the underlying continuum (see e.g. Schneider 1991), but the 
median value of 75 A is similar to that reported elsewhere for 
z ~ 4 QSOs, ~ 70 A (Fan et al 2001, Constantin et al 2002). 
However, 3 of our 18 radio-selected QSOs have rest-frame 
equivalent widths < 25 A : 0831+52 (< 10 A), 0918+06 (< 
10 A) and 1639+40 (~ 15 A), presumably because the lines 
are heavily self-absorbed. By comparison, only 5 out of 93 
optically-selected SDSS QSOs z > 3.6 (from Fan et al 1999a, 
Fan et al 2000, Fan et al 2001, Schneider et al 2001, Anderson 
et al 2001) exhibit such small Lya equivalent widths (but 
see Fan et al 1999b for an example of a z = 4.62 QSO with 
no emission lines at all). The difference between the two 
samples is weakly significant. Our sample of QSOs and that 
from SDSS have similar median redshift, so the difference 
cannot be ascribed to a difference in epoch (e.g. at z > 5.7, 
Fan et al 2003 found that 1 out of 3 QSOs has Lya equivalent 
width < 25 A). A difference might also arise due to the 
inverse correlation between emission-line equivalent width 
and optical luminosity known at low redshift, the Baldwin 
effect (Baldwin et al 1989). The optical luminosities of our 
QSOs are indeed ~ 2 mags higher than those from SDSS, 
but Constantin et al (2002) found the Baldwin effect to be 
weak or absent for QSOs with z > 4. Finally, the difference 
might be connected with the fact that our QSOs host a radio 
source, although previous studies (e.g. Corbin 1992) have 
found that emission-line strengths do not depend on whether 
a QSO is radio loud or quiet. If one combines our radio- 
selected sample (18 QSOs) with the radio-selected (but radio 
brighter) samples of Hook et al (2002; 13 QSOs z > 3.6 
with spectra) and of Snellen et al (2001; 4 QSOs z > 3.6, 
one already included in our sample), the fraction with Lya 
equivalent widths < 25 A reduces to 3 out of 34, i.e. not 
significantly different from SDSS. 

High-redshift QSOs are also characterised by a drop in 
the continuum across the Lya line, from red to blue, due to 
absorption by neutral hydrogen along the line of sight. The 
measured values of the opacity Da = (1 - fi / /i), where 
fi and f2 are the intensities in the continuum regions 1050 - 
1170 A and 1250 - 1350 A respectively (Oke & Korycansky 
1982), are given in column 16 of Table 1. The error in mea- 
sured D A is ~ 10%. The median D A for the 8 QSOs with 3.6 
< 2 < 4.0 is 0.50, while that for the 10 QSOs with 4.0 < z < 
4.4 is 0.62, consistent with the trend with redshift observed 
by Kennefick et al (1995b, their fig. 5) for optically-selected 
QSOs. 

4-1.3 Narrow absorption lines 

The QSOs 0831+52, 0918+06, 0941+51, 1303+00 and 1412- 
01 show metal-line absorption features at velocities relative 
to the QSO ranging 1200 to 7000 km s _1 . QSOs 0831+41, 
0918+06, 0941+51 and 1303+00 exhibit absorption at veloc- 
ity separations below 3000 km s _1 (i.e. probably associated 
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Figure 4. Distribution in SiAGHz and E mag of z > 3.6 radio- 
selected QSOs. • = FIRST QSOs (Table 1), S = Snellen et al 
(2001), H = Hook ct al (2002). Sec also Isaak et al (2002) for a 
list of 17 z > 4 optically- and X-ray-selected QSOs with known 
radio counterparts. 

with the QSO). 0831+41 and 1412-01 exhibit absorption 
systems at higher velocities. Metal absorption lines are not 
visible in our spectrum of 0747+27, but in a high-resolution 
spectrum obtained by Richards et al (2002), at least 14 in- 
dependent CIV absorption systems are detected. 

4-1-4 BAL QSOs 

One of the 18 z > 3.6 QSOs, 0831+41, is classified as a 
broad-absorption-line (BAL) QSO by Irwin et al (1998). Fan 
et al (2001) found 2 classical BALs in their sample of 39 high- 
z QSOs. The BAL fraction in both our sample and that of 
Fan et al is consistent with the fraction 10% found for low- 
redshift optically-selected QSOs (Weymann et al 1991). 

1624+37 (z = 3.4, i.e. not included in the sample of 
18 2 > 3.6 QSOs, see Table 1) is a new BAL QSO. The 
CIV absorption feature has FWHM ~ 9000 km s _1 and is 
detached by 20000 km s _1 from the peak of the emission. It 
has an unusually sharp high- velocity cutoff. 

4.2 Low-redshift objects 

The remaining candidates (Table 2) are a mixture of low- 
redshift QSOs, emission-line galaxies, absorption-line galax- 
ies, and stars. The QSOs and galaxies are mostly true iden- 
tifications, the stars are mostly misidentifications, as is im- 
plied by the distributions of radio-optical separations in Ta- 
ble 2 (although the higher values for stars might also be 
attributed to proper motion between the epochs of POSS-I 
and the the FIRST survey). The distribution of radio flux- 
densities of the sources apparently identified with stars is 
statistically indistinguishable from that of the FIRST cat- 
alogue as a whole. We obtained deep images at the JKT, 
R < 24, of two M stars lying close to Siaghz ~ 100 mjy 
sources (0741+33, E = 18.4, and 1353+34, E = 17.8 from 
Paper I), but a host galaxy is not detected in either case, 
so we cannot reject the hypothesis that these are true radio 
stars. 

A '?' classification in Table 2 indicates that there was 
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insufficient signal-to-noise to distinguish between galaxy and 
M-star spectra, but that the spectrum is not consistent with 
that of a 2 ~ 4 QSO (i.e. we expect there to be <~ 1 such 
QSO hidden amongst the '?' objects). 



5 CONCLUSIONS 

We identified 194 high-redshift QSO candidates which: 

(i) coincide with FIRST radio sources Siaghz > 1 mjy; 

(ii) are classified starlike on APM (POSS-I) and on 
POSS-II; 

(iii) have APM E < 18.8 and O-E > 3.0, or are invisible 
on the O plate. 

The sample covers an area of 7030 deg 2 . We have observed 
121 of the candidates and find 18 to be QSOs with z > 3.6, 
10 of these with z > 4 (Benn et al 2002, and this paper). We 
estimate that we have found 75% of the high-redshift QSOs 
present amongst the 194 candidates. The surface density of 
z > 4 QSOs with E < 18.8, Si iG Hz > 1 mjy, is 0.0031 
deg- 2 . 

This is currently the only well-defined sample of z ~ 4 
radio QSOs selected independently of radio spectral index. 
Vigotti et al (2003) use a subsample of 13 of these QSOs, 
with 2 > 3.8, M AS (1450 A) > -26.9 and log Piaghz(W 
Hz -1 ) > 25.7, to measure the space density of QSOs at z = 
4, and determine the change in space density between z = 
2 and 2 = 4. 

These QSOs are highly luminous in the optical (8 have 
M B < -28, q = 0.5, H = 50 kms^Mpc- 1 ). The SEDs are 
remarkably varied. They include: 3 QSOs with very weak 
Lya (0831+52, 0918+06, 1639+40), one with an unusually- 
high density of CIV absorption systems (0747+27, Richards 
et al 2002), one with a probable DLA (0941+51), the lumi- 
nous and much-studied 0831+52 (Irwin et al 1998), and a 
QSO with narrow Lya of high equivalent width (1309+57). 
1624+37 (outside the sample) is an unusual BAL QSO. 
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ins give: (1-2) optical right ascension and declination, (3) FIRST 1.4-GHz flux density (integrated, for Si. 4GHz > 10 m Jy> otherwise peak), (4,5) integrated Effelsberg 
and 10.6-GHz flux densities, (6) radio spectral index a, S v <x i> a , (7) logio radio luminosity (W Hz -1 ) (assuming H„ = 50 kms _1 Mpc _1 , q = 0.5, as used 
it in this paper, and assuming radio spectral index a = -0.3, as used by Vigotti et al 2003), (8) radio - optical position difference, (9-10) APM POSSI E magnitude 
3 colour, (11 - 12) JKT R-band and I-band magnitudes, (13-14) redshift and rms error on redshift, (15) approximate rest-frame Lya equivalent width (A), ± 30%, 
nuum decrement across Lya line (Section 4.1.2), (17) notes. The first section of the table lists the 8 z > 3.6 QSOs discovered amongst the 60 spectra reported here, 
ists the remaining 52 candidates, as well as 3 with redshifts from the literature, and 4 candidates classified as galaxies on the basis of JKT imaging.) The second 
the table lists the 10 z > 3.6 QSOs discovered amongst the 54 spectra reported by Benn et al (2002). The numbers 60, 3, 4 and 54 above sum to the total of 121 
5. The third section reports two z > 3 QSOs (with O — E < 3, not part of the sample) also discovered during the INT observations. 



